Electromagnetic Response of a Pinned Wigner Crystal 
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A microscopic model for analyzing the microwave absorption properties of a pinned, two-dimensional 
Wigner crystal in a strong perpendicular magnetic field is developed. The method focuses on excitations 
within the lowest Landau level, and corresponds to a quantum version of the harmonic approximation. 
For pure systems (no disorder), the method reproduces known results for the collective mode density of 
states of this system, and clearly identifies the origin of previously unexplained structure in this quantity. 
The application of the method to a simple diagonal disorder model uncovers a surprising result: a sharp 
(delta-function) response at zero temperature that is consistent with recent experiments. A simple spin 
lattice model is developed that reproduces the results of the quantum harmonic approximation, and 
shows that the sharp response is possible because the size scale Lc of patches moving together in the 
lowest frequency collective mode is extremely large compared to the sample size for physically relevant 
parameters. This result is found to be a direct repercussion of the long-range nature of the Coulomb 
interaction. Finally, the model is used to analyze different disorder potentials that may pin the Wigner 
crystal, and it is argued that interface disorder is likely to represent the dominant pinning source for 
the system. A simple model of the interface is shown to reproduce some of the experimental trends for 
the magnetic field dependence of the pinning resonance. 



PACS numbers: 73.40.Hm, 73.20.Mf, 75.40.Gb 
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I. INTRODUCTION 

It has long been appreciated that the groundstate 
of electrons in an otherwise structureless environment 
should be crystalline at low enough densities . Consid- 
erable effort has been focused on creating such a state in a 
two-dimensional electron gas (2DEG) as realized in semi- 
conductor heterojunction and quantum well systems [Q, 
although obtaining the appropriate limit of low electron 
and impurity densities has proven difficult. The intro- 
duction of a magnetic field perpendicular to the 2DEG 
improves this situation by raising the electron density 
at which crystallization is thought to occur. Recent ex- 
periments have identified the onset of strong insulating 
behavior in very high mobility systems, at magnetic fields 
such that the filling factor i' — N/N^, with N the number 
of electrons and the number of magnetic flux quanta 
through the system, is below ~ 1/5 for electrons |^,^ and 

1/3 for hole systems While much intriguing exper- 
imental data has accumulated, a definitive proof that the 
low filling factor insulating state of these two-dimensional 
systems is indeed an electron, or Wigner, crystal (WC) 
has remained elusive. 

One experimental fact for the low filling factor systems 
that is in clear agreement with expectations for a WC in- 
terpretation is that they are insulating. It is well appre- 
ciated by now, from analogous behavior in charge density 
wave (CDW) systems that an arbitrarily small disor- 
der potential should pin the WC at zero temperature so 
that there are no charge carriers that can flow in response 
to an arbitrarily small static electric fleld. In CDW sys- 
tems, the pinning potential also supplies a restoring force 
that induces a broad peak in the ac electromagnetic re- 



sponse at the "pinning frequency" , whose magnitude may 
be used to assess the correlation length of the CDW [^,0 . 
Early experiments on the magnetically-induced WC iden- 
tifled similar structure in the density response using var- 
ious techniques l^-[ll||. Theoretical analyses |]l2|-p^ of 
these experiments have largely focused on determining 
the correlation length of the WC from the experimental 
data, although no clear consensus on this quantity has 
yet been reached. 

Very recently, experiments on high quality hole 16 
and electron |17| systems at low filling factors have re- 
vealed structure in their microwave absorption prop- 
erties that are qualitatively different than what is ob- 
served in CDW systems. At the lowest temperatures 
and highest magnetic fields available, these systems ex- 
hibit sharp resonances at low frequencies, with quality 
factors Q = f / /S.f (/ here is the frequency of the reso- 
nance peak and A/ its width) as high as 30. The fre- 
quency of the peak increases with increasing magnetic 
field, and may saturate at a maximum for the largest 
fields ||l^. By contrast, most existing theoretical work 
predicts |l3,13 18 19] a broad {Q ^ 1) resonance, and 



a pinning frequency that decreases with magnetic fleld. 
The subject of this work is to understand some aspects 
of this puzzling experimental flnding. 

In what follows we will adopt a microscopic model of 
the magnetically induced two-dimensional Wigner crys- 
tal. The groundstate is assumed to be well described 
by a product of localized Gaussian wavepackets, so that 
exchange effects are ignored Our goal is to flnd 

the response of the system to a time dependent, spa- 
tially uniform electric fleld To compute the latter 
quantity we will employ a "quantum harmonic approxi- 
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mation" (QHA) |2^, a natural generalization of the clas- 
sical harmonic approximation, in which a finite number 
of angular momentum orbitals per site is retained in the 
Hilbert space of states for the electron system. This is of 
course sensible provided that the electrons remain close 
to their groundstate sites after being excited by the elec- 
tric field. A very strong magnetic field is assumed so 
that projection into the lowest Landau level is appropri- 
ate and all results reported here are for zero tem- 
perature. For an undisordered WC, the QHA reproduces 
known results found without resorting to a harmonic ap- 
proximation , and indeed offers an explanation of fea- 
tures found in the density response function that were not 
easily interpreted previously. 

Because exchange effects are negligible in this system 
, a time-dependent Hartree approximation turns out 
to be both convenient and accurate for computing the 
response functions of this system. For simplicity we con- 
sider only a diagonal disorder model; i.e., the disorder 
enters only in the single particle excitation energies and 
not in the interaction matrix elements among the elec- 
trons. In practice this means we have adopted a model 
in which the lattice is perfectly ordered, and disorder is 
introduced as a random on-site pinning potential |Q . As 
we will see, this assumption greatly simplifies the com- 
putation of matrix elements entering the QHA. Using a 
perturbative approach described below, for weak disorder 
it is possible to demonstrate that our qualitative results 
are insensitive to the diagonal disorder assumption. 

For all the disorder models studied in this work, the 
results of the QHA turn out to be qualitatively the same: 
the response functions exhibit an extremely sharp re- 
sponse at the lowest collective mode frequency (provided 
that the electron-electron interaction is unscreened; see 
below.) This is the central result of this work. Two phys- 
ical models of pinning are analyzed in detail in this study. 
Charged impurities which may be found in the spacer 
layer between donors and the electrons have been argued 
p6| to represent the strongest source of pinning for the 
WC in dc non-linear I — V measurements; i.e., these set 
the energy scale necessary to fully dislodge the electrons 
in a non-vanishing static electric field so that a current 
may flow. From the CDW point of view 0, this model 
corresponds to a strongly pinned system. However, we 
will see that the strong pinning approach grossly over- 
estimates the pinning frequency, so that the experimen- 
tally measured pinning frequencies cannot be explained 
by this source of disorder The ultimate reason for 

the discrepency between the pinning frequency in the mi- 
croscopic model and the CDW result for this is that a 
magnetically-induced WC is not a CDW. Our results in- 
dicate that the collective mode spectrum associated with 
strongly pinned centers is far more like that of a crys- 
tal with vacancies pinned at charged impurity sites than 
that of an elastic medium tied down at random sites. 
We then focus on a model of interface disorder that 



we believe is likely to be the dominant pinning mech- 
anism in electromagnetic absorption for these systems. 
Heterostructure and quantum well interfaces are believed 
to have structure at length scales of several tens of 
angstroms . An interpretation of this is that in semi- 
conductor hcterostructures the interface between differ- 
ent materials (typically GaAs and AlAs) may only be 
defined to within one or two lattice constants (~ 5 A). 
This is often modeled as an interface with pits and/or 
islands of typical size scale < 100 A Electrons re- 

siding in large pits or regions with an unusually large 
number of small pits have an enhanced probability of 
lying slightly closer to the donor layer and so may be 
bound by them. An important aspect of the physics in 
this model is that in the lowest Landau level interfaces 
with pits whose size scale are smaller than the magnetic 
length ^0 = {hc/eBy/^ have pinning potentials that in- 
crease with decreasing Iq. This leads to an increasing 
pinning frequency with magnetic field, as seen in experi- 
ment. The ~1 GHz magnitude of the measured pinning 
frequency may be explained with very reasonable param- 
eters describing the interface at which the 2DEG resides, 
as described below. 

A typical result for a single disorder realization, whose 
precise form is described below, is shown in Fig. 1. The 
system includes 1024 electrons, and is assumed here (as 
throughout this work) to obey periodic boundary condi- 
tions. Because of the finite number of degrees of freedom, 
the density response function consists of a series of delta 
functions. Note that the scale of the figure is logarith- 
mic, so that a linear plot of power absorbed as a func- 
tion of frequency shows a single delta-function response 
at the lowest collective mode frequency, as illustrated in 
the inset. The height of this peak shows no sign of de- 
creasing with increasing system size, nor do any other 
collective modes develop a noticeable oscillator strength 
in power absorption. For a model in which the typical 
on-site pinning potential is small, we find the pinning fre- 
quency (i.e., the position of the delta- function response 
in Fig. 1) to a good approximation is given by vq, the 
energy required to excite an electron from the m — 
to the TO = 1 angular momentum state in the absence 
of other electrons, averaged over all the electrons. This 
value is is essentially the same as the well-known result 
for weakly pinned CDW's 0|. Typically, the exact pin- 
ning frequency found in the QHA falls somewhat below 
this estimate, by an amount that is specific to the precise 
disorder realization. This correction to the weak pinning 
result may be estimated using a perturbative approach 
described below. 

Although the sharpness of the absorption peak seems 
to be in agreement with recent experiments, one needs 
some analytic method to demonstrate that in the ther- 
modynamic limit the peak does not broaden, particularly 
since this result is so different than prior expectations 
Towards this end we develop a (pseudo)spin lat- 
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tice model in which the z component of the spin at a 
given site represents the angular momentum state of an 
electron. A convenient perturbation theory for the sys- 
tem may be developed around a uniformly pinned state 
(i.e., one in which the pinning potential for every site is 
1^0.) The perturbing parameter is then AU{R), with {R} 
representing the sites around which the electrons are lo- 
calized and U{R) = vq + AU{R) the energy required to 
excite an electron from the m = to the m = 1 state. For 
the uniformly pinned system, it is natural for the power 
absorbed to have a (5-function response, since the collec- 
tive modes have wavevector g as a good quantum number. 
The correction to the weight W of this (5-function due to 
disorder may be computed in perturbation theory, and a 
depletion D may defined such that W oc 1 — D. To low- 
est non-trivial order in perturbation theory, the largest 
contribution to D is found to have the form: 



D oc 



\AUiq)\'\u+\' 
\El~ |2 



(1) 



Here q represents wavevector, AU{q) the Fourier trans- 
form of the perturbation, and E^ the collective mode 
dispersion for the uniformly pinned system. is a 

weighting function whose precise form is given below; in 
the limit q 0, \ut\ 1. Provided that the disorder 
potential does not have fluctuations on arbitrarily long 
length scales - i.e., that arbitrarily large patches of un- 
pinned or strongly pinned electrons are rare - then one 
expects \AU{q)\'^ oc q^ for small q. The dispersion rela- 
tion for the uniformly pinned system we show below has 
the form 



EP^v^ + 2Tie'pnllqlK + 0{q^) 



(2) 



for small q, where po is the electron density and k, is the 
dielectric constant of the host semiconductor. The lin- 
ear dispersion with q turns out to be a direct result of 
the long-range nature of the Coulomb interaction, and is 
not present in a model with short-range (e.g., screened) 
electron-electron interactions. Plugging Eq. || into Eq. 
|l| demonstrates that for a given disorder realization a fi- 
nite depletion results that will be small if the pinning 
is not too far from uniform. This directly demonstrates 
that it is the long-range nature of the electron-electron 
interaction that is responsible for the sharp response in 
the system. A system with short-range interactions has 
a collective mode spectrum dispersing quadratically with 
q away from vq, and as seen for Eq. |l|, this leads to a 
divergence in D, signaling that the response is not sharp 
in this case. Fig. 2 illustrates the response function 
for a model in which the electron-electron interaction is 
screened. As may be seen, even for a relatively small 
number of particles {N = 225), the response has been 
significantly broadened, as is typically found in CDW 
systems [0,0,01 . Roughly speaking, this says that an 



appropriate analogy for the randomly pinned WC in a 
strong magnetic field is a hard object that does not de- 
form much when vibrating on randomly placed springs, 
so that a well-defined periodic response is to be expected 
in spite of the randomness of the pinning. 

One important caveat to this result must be noted. 
The assumption that |AC/(g)p oc q^ is actually not true 
for the most common forms of disorder studied. Specifi- 
cally, white noise pinning models usually have equal fiuc- 
tuation strengths at all length scales, so that AU does 
not vanish as g — > 0. In our expression Eq. ^ this leads 
to a formal logarithmic divergence in the thermodynamic 
limit and an expected broadening of the absorption peak. 
We note however that the disorder models used in the 
QHA have a white-noise character, yet we have not been 
able to detect a finite width in the resonance for any 
disorder realization. Presumably this indicates that the 
(5-function response is a finite system size effect, and one 
needs to estimate a length scale above which broadening 
in the response function should become apparent. This 
may be accomplished by noting that the integral in Eq. 
|l| has an infrared cutoff 27r/L where L is the system size. 
The value of L for which 13 « 1, which we call Lc, sets 
the size scale above which broadening will be significant. 
Solving for Lc, one finds 



^ exp{ 



\AU[ 



where is a length scale above which Eq. || is accurate 
[£ Ri lOflo, with ao the WC lattice constant], and AU = 
liniq^o AU{q)K.lo/e^. 

The physical interpretation of Lc is that it is the length 
scale for coherent motion of patches of electrons in the 
lowest collective mode, since for L ^ Lc the perturba- 
tion theory is valid and the system responds much like as 
a single oscillator to the excitation field. A conservative 
estimate of Lc for the experiments of Refs. shows 
that in practice it is huge, Lc > lO^'^ao &tv = 0.22, where 
flo is the lattice constant of the WC. Thus, remarkably, 
in spite of the formal divergence in the thermodynamic 
limit, the perturbation theory is in fact controlled and 
valid for samples with physically relevant dimensions, 
where L ~ IO^oq ^ Lc- Thus, our interpretation of 
the WC response as one of an undeformable oscillator is 
indeed appropriate. 

This article is organized as follows. In Section II, we 
develop the QHA used to compute the electromagnetic 
response of the pinned WC. Section III develops the pseu- 
dospin approach to the collective modes and shows how 
a perturbation theory may be developed around the uni- 
formly pinned state, and the calculation of the depletion 
D is discussed. Section IV discusses some details of the 
interface pinning model. Readers interested in our re- 
sults and not the details of the calculations may wish to 
proceed directly to Section V, which discusses results for 



3 



different pinning models, including charged impurities in 
the electron layer and the interface pinning model. Sec- 
tion VI discusses the relationship of this work with results 
of other calculations, as well as several unresolved ques- 
tions regarding the experimental results. We conclude 
with a summary in Section VII. 



II. QUANTUM HARMONIC APPROXIMATION 

In this section we will develop a quantum mechanical 
generalization of the harmonic approximation, in which 
the electrons are treated as distinquishable, so that the 
groundstate and low-lying excited states may be repre- 
sented as linear combinations of (unsymmetrized) prod- 
ucts of single particle states. The single particle states 
we will use consist of angular momentum states, and we 
will find that it is an excellent approximation to retain 
only the lowest lying ones when computing low-energy 
properties of the system. Once we have truncated the 
Hilbert space in this way, it becomes possible to numer- 
ically compute response functions for reasonably large 
systems, which one may compare with experiment. 

The groundstate of two-dimensional charged particles 
in a very strong magnetic field is believed to be accu- 
rately represented by a collection of distinquishable par- 
ticles in Gaussian orbitals of the form 
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-\fi-R.\^/4ll+^z-(RiXri)/ll 



In this state, the electron with coordinate rt is maximally 
localized within the lowest Landau level. Since the ki- 
netic energy of a collection of electrons in such orbitals is 
already as small as possible, the groundstate presumably 
will be found by minimizing the total potential energy 
with respect to the parameters Ri. In the absence of dis- 
order these should be chosen to lie on a perfect lattice, so 
that the charge distribution in the limit of infinite mag- 
netic field {Iq ~^ 0) approaches that of the groundstate 
for a distribution of classical point particles. Disorder 
changes the optimal positions for the electrons |^^, but 
provided the disorder is not too strong, the site centers Ri 
will presumably not be close to one another on the scale 
of the magnetic length Iq. In this situation, wavefunction 
overlaps among the individual electrons are negligible, so 
that a product of the single particle wavepackets (j)Q j^.i^i) 
is an accurate representation of the groundstate, despite 
the fact that it is not explicitly antisymmetric pc| ]. 

Since we are interested in the electromagnetic response 
of the system at low frequencies, it will be necessary to 
compute the low-lying excited states of the system. In 
a classical analysis, this is typically done by assuming 
the displacements of the electrons from their groundstate 
positions are small, so that the energy may be expanded 
to second order in displacements, allowing the normal 



modes of the system to be found in a relatively simple 
fashion. If we assume the electrons are localized within 
a magnetic length of their groundstate positions, then 
the analog of the small displacements approximation is 
to allow individual electrons to be excited into higher 
angular momentum states localized about each lattice 
site. Provided the angular momentum is not too large, 
the charge densities associated with these states for a 
given site will not overlap significantly with those of its 
neighbors, so that exchange effects may continue to be 
ignored. Thus we consider a set of states for each site 
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,1/2 



27r/g2™m! 

,-|r-i?.|V4i= + J5-(i?,xr)/i^ 



(3) 



where z ~ x + iy is the electron position in complex nota- 
tion, and Zi — Rf + iR^. If we define creation operators 
for these single particle states al^^, then the Hamiltonian 
for the system may approximately be written in the form 



H 



-It. e 

i^j minirrijnj 



minim .in 



Here, represents the interaction of an isolated elec- 
tron with a disorder potential, and the interaction matrix 
element is given by 



^ mim2m3m4 



/ j2 j2 { ^ \'^{ ^ -1 (mi+m2+m3+m4)/2 







L2/2J 

ZtQ 



xz*'"iz^"^Z2*"^z™*e-'^i/2'o-'-2/2i^„(^^ „ + R^j), (4) 

where for most cases we will study, v{r) is the Coulomb 
interaction e^/^r, with k the dielectric constant of the 
host material for the electron layer. 

To simplify our calculations, we will make certain as- 
sumptions about the form of the disorder entering H . 
Firstly, we assume that although disorder certainly will 
cause the orbit centers Ri to vary from the positions of a 
perfect lattice, this variation does not affect the qualita- 
tive features of the absorption spectrum. Thus in prac- 
tice our RiS take on the values of perfect triangular lat- 
tice positions. It should be noted that this choice of the 
RiS is not required to carry out the QHA, but by adopt- 
ing it the numerical computation of U^^.^.^.^^. is greatly 
simplified. Some details on how this is done in practice 
are discussed in Appendix A. The second simplification 
we introduce is to assume that the pinning potential at 
the individual sites is circularly symmetric, so that in the 
groundstate the electrons occupy the 4>m=o,i orbitals and 
do not admix higher angular momentum states. Again, 
the QHA may be developed without this assumption. 
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However, this simpHfication has the great advantage of 
allowing an analytic specification of the groundstate. In 
the absence of this assumption one would need to find 
the groundstate orbital occupations numerically, which 
presumably could be computed with sufficient accuracy 
using a static Hartree approximation. While such a cal- 
culation is clearly feasible, we believe that allowing as- 
symmetric forms for will have little quantitative ef- 
fect, provided the average of this quantity over the sites 
restores the circular symmetry. Finally, to take full ad- 
vantage of the symmetries of the system, we will impose 
periodic boundary conditions. 

A second caviat related to this is that in principle one 
must retain fewer than seven orbitals on each site to 
have each electron lie purely in the m = orbital in 
the groundstate. Even in the absence of disorder, the 
six-fold symmetry of the lattice will in principle admix 
in states with angular momenta equal to integral multi- 
ples of six. Although this effect is extremely small [ po[ , 
it can in principle lead to weak instabilities in the re- 
sponse function computed below when too many orbitals 
are included per site. As we will see, all the low energy 
excitations (i.e., phonons) of the system are accurately 
captured when only two orbitals per site are retained, so 
that in practice this does not arise as a problem. 

With these simplifications, the Hamiltonian takes the 
form 



i m 



(5) 



i^j minim jUj 



The zero of energy may always be chosen such that 
ei{m = 0) = 0, and we expect ei{m > 0) > 0. The 
choice of {ei(m)} defines the specific disorder model we 
are studying. We note that the greatest power of this 
method is that it is well-suited to disorder potentials that 
vary on length scales of the order of Zq and smaller. The 
two pinning mechanisms we study in detail - charged 
impurities and interface pinning - both fall into this cat- 
egory. 

The quantity we will calculate, from which either 
power absorption or frequency dependent conductivity 
may be computed, is the response function 

(6) 

where here t is imaginary time and amiij), aj„j(''') are 
the usual time-dependent Heisenberg representations of 
the annihilation and creation operators, and the time 
ordering operator . The brackets < • • • > represents a 
thermal average. The Fourier transform of this function 
with respect to imaginary time has poles at the collective 
mode frequencies. To generate a closed formula for x, we 



employ an equation of motion method similar to that 
used in Ref. [Q. The time derivative of x satisfies the 
equation 



rrL2'm3m4 



dT 



= - < [Pmim2(*;'^).Pm3m4 0';0)] > K^) 

(7) 



< Pm^mA3\^) > 



dT 

where Pmim2{^T) = a^Jrl^^iT)am2^{'r) ■ The time deriva- 
tive of p may be computed from the Heisenberg equation 
of motion 



Note that the formal inclusion of the chemical potential 
/i is necessary because of the use of the finite temperature 
formalism; this means our formalism allows in principle 
for us to treat the situation in which there is more than 
one electron per site. In what follows, /i will be chosen 
such that there is one electron per site in the groundstate, 
and we will see in any case that any formal dependence on 
/i drops out of our equations for the response functions, 
so that it does not need to be explicitly computed. The 
relevant commutators in the equation of motion may be 
worked out using {ami,a}^j\ = SmnSij from which 
one may show 



/^m 17714 m2m3 Pm^^m^ {j'^^'inim^)^!^ ■ 

(8) 

Combining Eqs. ||, 0, and|8[ the equation of motion takes 
the form 

dT 

Xmim2m3m4 



E 

nin2"3"4 I 



^ ^ ^rti n^n^m < [Pnim2 (*i ''')Pn3n4 (^5 ■'")'5n2mi 



> 



(j;0)>. (9) 

Here C/;;i,„3„, = U'^^^,^„.^„^{1 - Su). To obtain a closed 
form for the equation of motion, we utilize a Hartree de- 
composition of the last term in Eq. ^. This is very much 
in the spirit of Ref. [ p4[ ; however, because we are ignoring 
overlap among single particle states located at different 
sites, it is unnecessary (and in fact would be inappropri- 
ate) to include exchange terms in the decomposition. We 
thus make the substitution 

(Pnim2 (*i T)Pn3ni (^i '^)Sn2mi Pmin2 (*i '^^Pn^rii 'T)Snim2'} 
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~t~(*^ Pnim2 (^) ^n2mi ^ Pmxn^i.'^^ ^ ^ni m2 ) Pn3n4 ; • 

Substituting the above decomposition and Fourier trans- 
forming Eq. ^ with respect to imaginary time, we arrive 
at the time-dependent Hartree approximation for the re- 
sponse function: 



"("^ Pmim4(^) ^ ^m2rii3 Pm3m2 (^) ^ ^77117714) 



+ (£i(rni) - £i(m2)) 



Xmir 



ni 7712 



77371477737774 



(0'; j^^n) 



(10) 



Once we have solved Eq. |l^, in principle we may compute 
any response function we like. However, since ultimately 
we are interested in the conductivity or power absorp- 
tion of the system for a spatially uniform electric field, 
it is convenient to formulate the equation for a response 
function that is less cumbersome in terms of indices, but 
which may still be used to compute the physical quan- 
tities of interest. Towards this end we just consider the 
center of mass response of the system. An operator cor- 
responding to the displacement of the center of mass may 
be written in the form 



UCM 



^ 'Y < mi,i\r- Ri\ni2,i > al^^^am2,t, 
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where \m2,i > is a ket-vector representation of 4'm,i{'r)- 
The perturbation due to a time-dependent, spatially uni- 
form electric field (e.g., microwaves) in terms of ucm is 
—cNEq ■ Ucm, where N is the total number of electrons. 
A quantity whose response to this perturbation is conve- 
nient to study is the center of mass displacement itself, 
so that the response function we will actually focus on is 



(12) 



The conductivity and hence power absorption may 
be written in terms of this quantity as follows. The 
time Fourier transform of the spatially averaged cur- 
rent density in response to an applied external field 
Re-Eoe'"* is < j{uj) >— iepoOJUCMioj), with po the 
sheet density of the electrons. In linear response theory, 
this takes the form (at zero temperature) < j{uj) >= 

—ie^ PqojN I^Xcfh-ii^ + *i5)^'o: where in the usual way 
pl| we have Fourier transformed the response function 
with respect to imaginary time, and made the replace- 
ment iujn io + i5 to take the zero temperature limit. 
The dot product of this quantity with the total elec- 
tric field is proportional to the power absorbed. A mi- 
nor complication is that one must include the screen- 
ing field generated by the displacements of the elec- 
trons. Since we are interested in the bulk current through 
the system, we assume that the induced electric field 
may be replaced by its spatial average Q. Using 
a dipole approximation for the electric field generated 
by the motion of the electrons, it is easily shown that 
Eind{t) = ^aucM{t), with a = J2i^o^/^i- Using 
Etot ^ Eq + Eind, one finds for the power absorption 
< i'(^) > -Etot oc Y.a.,0^tot^a(3{uj)E^^^, with & the con- 
ductivity matrix given by 

a{uj) = -ie^pQU}NxcM{^ + - 2aNxcM{^ + 

where here we are regarding a and xcm as 2 x 2 ma- 
trices. As has been pointed out before Q, the induced 
electric field shifts the frequency of the peaks in a{ijj) 
from where they are found in xcm(w -|- id). However, 
since the width of the peak we will find in xcm is re- 
markably small, this shift may be neglected for our pur- 
poses. Thus in this study we will focus on xcm ^nd its 
frequency dependence. 

To compute the center of mass response, we define an 
intermediate response function 



Xr. 



in terms of which 



I mim2 



Combining the definition of Xmima (*'''") "^ith Eq- |l^, 
one derives the equation of motion 



where a, /3 = a;, y. 
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< Pmim(i) >< TO2|r'^|m >] 



|lj may be written schematically in the form 

i,mi,m2),{l,ni,n2) 

ni ,n2 ,1 



+ [£j(mi) - e,(m2)] Xmims *^») 



^n2''7ii Prnin2 

{^)> 



Note that the matrix elements < 



(13) 

r^\mi > are the same 
as those appearing in Eq. 11, evaluated for Ri — 0. 

Eq. |l^ may be solved if one knows the groundstate 
densities < Pmimiii) >■ It is here that our assump- 
tion that all the electrons reside fully in the lowest an- 
gular momentum orbitals for their sites in the ground- 
state becomes useful: the form of the densities is simply 
< Pmim2{i) >= <5miO'^m20- Furthermore, as shown in Ap- 



ex 5„ 



pendix A, it can be proven that U'^l 
that Eq. O is considerably simplified, taking the form 



I il 







< 0\r'^\mi > S^,_o- < m2\r'^\0 > 5m,o 



-|-[ej(mi) -£,(m2)]Xmi,„2(«;«^n) 



u: 



n^Ui I 



(14) 



where ei{m) — Si{m)+^i U'^^qq. It is also easy to show 
that Xmim2 = if either mi and m.2 are both zero, or if 
both are non-zero. This is the equation we actually work 
with, and it is solved in a way closely analogous to the 
method developed in Ref. HJ]. We regard Xmim2(0 a 
vector whose components are labeled by TOi,m2,L Eq. 



-^^(imim2),(i"in2)] Xnin^ (^i 



where A£i(mi,m2) = ei(mi) — Si{'ni2) is the energy re- 
quired to excite an electron on site i from the mi orbital 
into the m2 orbital when all the other electrons are static. 



U, 



D 

(imim2),('ni"2) 



HI ST 

0miTi3n4"m2:0 



u: 



HI r 
m20n3n4"'ni,0 



is the change in this energy difference due to the 
fact that all the electrons are in fact dynamic, and 

Xmrm2(*;*^n) 0|r'^|TOl > 5,„20- < m2\r^\Q > S„i^a 

is the response function in the absence of both electron- 
electron interactions and the disorder potential. The ma- 
trix M may be diagonalized numerically; it will have real 
eigenvalues LOj (which may be shown to come in pairs of 
equal magnitudes but opposite signs) and eigenvectors 
{mim2l). It can be shown that when regarded as a 
matrix in the indices j and {mi'm2l), V has an inverse 
. Denoting ujj as the diagonal matrix of eigenvalues 
of A/, the solution to Eq. 



14 



may be written schemati- 



cally as 



From this form, it may be seen that the poles of x, and 
hence the energies of collective modes of the system, are 
given by the set of eigenvalues LHj . 

As a first example, as well as a check on whether this 
method works, we consider the situation of a WC in an 
undisordered environment. One only needs then to set 
Si{m) = for all i, m in Eq. and proceed as de- 
scribed above. A histogram of the collective mode ener- 
gies, which represents the density of states for the sys- 
tem, is shown in Fig. 3. In this calculation, the number 
of electrons was N = 225, the filling factor was set to 
v = 0.2, and 5 states per site were retained. As may 
be seen, a broad peak near the origin is accompanied by 
three larger, well-separated peaks at higher energy. The 
number of such sharp peaks appearing in the density of 
states we find in general to be equal to the number of 
orbitals per site above the groundstate retained in the 
calculation. 

The meaning of the sharp peaks may be interpreted if 
one notes that their positions are very close to the values 
of Aei(m, 0) used in the equations of motion. This means 
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that the collective modes are only very weakly affected 
by Ud , so that on may understand the modes as local ex- 
citations in which an electron is excited to a high angular 
momentum. Thus, these peaks may properly be under- 
stood as an analog of Wannier excitons in tight-binding 
models of electron systems. Interestingly, these peaks 
have been previously noted in an approach that does not 
use a harmonic approximation [Q. However, in that ap- 
proach, the nature of these excitations was unclear; the 
introduction of the harmonic approximation allows us to 
understand why there is structure in the density of states 
at these energies. A careful comparison of the peak po- 
sitions for V = 0.25 found in Ref. (see Figs. 2 and 
3 of that work) with results from the method described 
here shows that the energies of these peaks are nearly 
identical for the lowest energy peaks. 

Fig. 4 illustrates the density of states for a calculation 
with N ~ 529, v — 0.2, and just two states per sites 
retained. As may be seen, the density of states is pre- 
cisely what one expects for a phonon density of states in 
a magnetic field: the leading edge of the phonon density 
of states rises sharply, consistent with the D{u!) oc ui^^^ 
expected for the low-frequency behavior of the density 
of states that results from the uj{k) cx A:'^/^ dispersion 
of the WC in a magnetic field . The peak contains 
a strong cusp structure consistent with a van Hove sin- 
gularity, and the width of the peak is found to decrease 
with increasing magnetic field, consistent with the l/B 
dropoff expected for the phonon bandwidth . 

The quantitative agreement of these results with pre- 
vious calculations give a strong indication that the QHA 
works, and we learn from the above calculations that all 
the low-energy modes (i.e., phonons) of the WC can be 
captured by retaining just two states per site. This con- 
siderable simplification allows one to treat fairly large 
numbers of electrons for the pinned WC; our largest cal- 
culations contain = 1024 electrons. From this point 
onward we will adopt the two-state approximation, and 
denote the energy difference between the two states on a 
site as Aei{m — 1,0) as Ae^ 

Having established that the QHA gives sensible results 
for cases where the correct answer is known, we can now 
use it to investigate the case of the pinned WC. We defer 
detailed discussions of the results to Section 0; however, 
they can be summarized extremely briefly. Essentially all 
of our center of mass response functions are dominated 
by a single collective mode, the lowest excited state aris- 
ing for any given disorder realization (e.g. Fig. 1). For 
the weak disorder models [i.e., Asi <^ ujb, where is 
the bandwidth of the phonon density of states (Fig. 4)], 
the position of the pinning frequency turns out to be ex- 
tremely close to V{), the average pinning energy per site. 
We find this form of the center of mass response for all 
our disorder models and system sizes studied, and there 
is no noticeably trend for the strength of this one mode 
to decrease with increasing system size within the sizes 



one may study using this model. One is led to conclude 
that the electromagnetic response of the WC in a strong 
magnetic field at zero temperature is extremely sharp, 
consistent with experiment but not with previous theo- 
retical expectations. As explained in the introduction, 
this turns out to be a consequence of the long-range na- 
ture of the Coulomb interaction. To see precisely how 
this result arises, as well as to assess whether it will sur- 
vive in the thermodynamic limit, we need to develop an 
approach from which one may learn what happens to this 
mode for very large systems. This is the subject of the 
next section. 

III. PSEUDOSPIN REPRESENTATION 

The results of the QHA indicate that generically there 
is a sharp resonance in the low frequency response of the 
magnetically induced WC for the finite size systems one 
is able to handle with that method. While there is no 
indication within those calculations that the resonance 
either broadens or weakens with increasing system size, 
one cannot rule out based upon them the possibility that 
a broadening does develop at very large system sizes. 
Thus one would like to develop models that are analyti- 
cally tractable wherein the thermodynamic limit may be 
taken. 

A. Quantum Spin Model 

In this subsection, we will develop a mapping of the 
electron system onto an effective spin lattice system, and 
obtain couplings between these spins using an expansion 
in Zo/ftQj with oq the nearest neighbor distance. One 
important advantage of this formulation is that the re- 
sulting couplings are far more tractable than the matrix 
elements found in the QHA, allowing us to make con- 
siderable progress analytically. The spin waves of the 
model are the analog of the phonons for the WC, and it 
will be shown that the resulting dispersion of these spin 
waves is identical to the classical phonon spectrum of the 
WC, demonstrating that the mapping produces sensible 
results. 

We begin with the observation from the QHA analy- 
sis that essentially all the low energy excitations of the 
magnetically-induced WC may be obtained in a model 
that retains only two states per lattice site. This moti- 
vates an approach in which one may wish to consider the 
on-site degrees of freedom of each electron as an effective 
pseudospin, with the m = angular momentum state 
representing a "spin up" state, and the m = 1 a "spin 
down" state. (Note that we will assume the real spins 
of the electrons are polarized by the magnetic field and 
do not represent a low-energy degree of freedom for the 
system.) The system is then mapped onto a quantum 
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magnet, whose interactions we will see are of a magnetic 
dipole form, and whose spin waves will have the same 
dispersion relation as the phonon spectrum of the under- 
lying classical electron degrees of freedom. Formally, the 
mapping is given by: 



25? 



(15) 



To write an effective Hamiltonian for these pseudospin 
degrees of freedom, it is convenient to employ a multipole 
expansion of the electron density, which becomes quan- 
titatively accurate in the limit lo/ao ^ 1 (i.e., v <C 1), 
where oq is the nearest neighbor distance. Neglecting 
exchange effects, the interaction energy is formally 



2 ^ 



,2 ,2 p{ri+Rt)p{r2+Rj) 

d rid r2 —. —■ 

\ri - r2 



(16) 



where a,/3 = x,y, J\ 



a, (3 _ 



hi = - iw' + = 0) - ^ii^n = 1) is an ef- 

fective magnetic field, which may be non-uniform due to 
the pinning potential. Note that we have dropped an 
irrelevant constant from the Hamiltonian. 

The low-energy collective modes of this model are spin 
waves, and in the absence of a pinning potential their dis- 
persion relation should be identical to that of the phonon 
modes of the underlying WC from which the pseudospin 
model was derived. We thus begin our analysis by consid- 
ering the case ei{m) = 0. To derive the spin- wave spec- 
trum, we rewrite the spin operators in terms of bosonic 
degrees of freedom p& with the approximate mapping: 



and 



-b!b, 



(19) 



Writing p{r + Ri) = pi{r)^ the densities may be written 
in terms of the pseudospin operators. 



+0S(f)0i(f)5- + ^Ur)Mr)S+, (17) 

where 0o(i) ^-re the m = 0(1) orbitals localized around 
the site of interest (Eq. |[) , and we have suppressed the 
exphcit site labels Substituting Eq. |l^ into Eq. 

P^ , expanding the integrand to second order in ri(^2)/Rij 
where Rij = \Ri — Rj\, and performing the integra- 
tions, one arrives at the formula U ~ ^ X^i^j i with 



' Ri 



2i?? 



2eH 



+—Jl[sls'j-3{SlRi,)iSlRi,)/R:, 



i?3 



where S [ 



j.^.-^ = {S^^-yS^^f^-yQ). Note that the effective in- 
teraction between spins in Uij is of an XY dipole form, 
which is not surprising given the fact that a small r ex- 
pansion for the electron-electron interaction leads to an 
electric dipole interaction. The third term in Uij takes 
the form of an effective magnetic field that tends to orient 
the spins in the +z direction; this reflects the fact that 
the <^o state is the groundstate of an electron in a given 
site if all the other electrons are fixed in their (/ig states. 
The Hamiltonian for the effective spin system thus may 
be written in the form 



a,f3 i^j 



(18) 



This mapping is an approximate form of the Holstein- 
Primakoff transformation | |3^ , expanded for the situa- 
tion < blbi >^ 1/2, where < • • • > is an expectation 
value for any of the low energy states that are of inter- 
est in the zero temperature response function. A simple 
way to see that < blbi > is small is to write it directly 
in terms of the underlying electron creation and anni- 
hilation operators: foj = a\^aQi, bi = a\^aii. One then 
finds b\] = a^iOiQi — a\^aii. If the mapping were exact, 
this commutator would be unity. However, to the extent 
that the groundstate is well-approximated by electrons in 



Gaussian orbitals, < aJjOoi >= 1 and < a\^aii >= in 



the groundstate. Furthermore, the lowest excited states 
- the spin waves in the bosonic language - are described 
in terms of the electronic degrees of freedom by a single 
particle excited out of the to = state into the to = 1 
state, averaged with an appropriate phase factor over all 
the sites. Thus, the expectation value of for the 

low-lying states is the same as that of the groundstate, up 
to corrections that vanish in the thermodynamic limit. It 
should be noted that when quantum (or, at finite tem- 
perature, thermal) fluctuations are important, then one 
must work with the exact mapping between spin and bo- 
son operators |3^], and higher order terms in bjbi should 
be retained, introducing spin wave interactions. Such 
terms are the analog of anharmonic terms in the Hamil- 
tonian for the underlying WC degrees of freedom, and 
we expect them to be very small at zero temperature for 
large magnetic fields (small , where placing the electron 
in Gaussian orbitals at specified locations that minimize 
the potential energy is thought to be an excellent approx- 
imation for the groundstate pO|] . From this point onward 
we will ignore spin-wave interactions, and our goal will be 
to understand why disorder does not broaden the electro- 
magnetic response associated with the quadratic Hamil- 
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tonian below. 

In terms of the bosonic operators, the Hamiltonian we 
thus will consider takes the form 

i 




Here, Ui = —hi, and = (i?f^ — iR\-)l Rij is a complex 
representation of the vector direction separating sites Ri 
and Rj. In the absence of pinning [ei{m) — 0], H is 
diagonalized in two steps. First, a canonical transfor- 
mation from real to momentum space separates out the 
independent modes: 

k 

-~Y.[G{k)bj:b_~, + G*{k)bl.b\], 

k 

where 6,- = F{k) = 1^^.41(1- 

e-*^0, and G{k) = 3 ^(n^p2g»fc.i?,^ The 

diagonalization of the Hamiltonian is then completed 
with a Bogoliubov transformation of the form 7^: = 
u^b^ + vgb^^^, with — l^gP = 1 guaranteeing that 
[Tfc' "^ll ~ "^^^'^ boson Hamiltonian then may be written 
as H ^ Fipl\lk '^^ ^^'^ chooses 

= cosh%e*'^£/2 
vj: = sinh^e-^-^E/^ 

E^^ = ^ F{kY - A\G{k)V (21) 

with G(k) = e*'^s|G(fc)|, and tanh26'j: = -2\G{k)\l F(k). 

thus represents the spin-wave dispersion for this 
model, and we need to confirm that it has the same form 
as the underlying electron degrees of freedom. For the 
WC in two dimensions in the absence of a magnetic field, 
the phonon dispersion is found by solving the eigenvalue 
equation j38| 

Y,C^p{k)ep^uj''{k)ep, (22) 


where 

J J J 

(23) 



and m* is the effective mass of the electrons. For every 
value of fc, Eq. ^ has two eigenvalues, corresponding; to a 
longitudinal mode Wi(fc) and a transverse mode ujt{k). In 
the presence of a magnetic field, the equations of motion 
for the electrons mix these two modes. One then obtains 
two different normal modes, one dispersing from the cy- 
clotron frequency ujc — eB/m*c, and the other having the 
form in a the strong field limit uj{k) = uji{k)ujt{k)/ujc- 
By solving Eq. |2^ and using Eq. one may show with 
some algebra that uii(k)ujt{k)/(jJc = E^. Thus the spin 
waves of our pseusospin model faithfully reproduce the 
exact phonon spectrum for the WC in a strong magnetic 
field, and we see that our model quantitatively captures 
the low energy dynamics of the WC. 

B. Uniformly Pinned Wigner Crystal 

In this subsection we will analyze the response function 
of a WC in which each site has precisely the same pinning 
potential. We will find that the power absorption from a 
spatially uniform, time- varying electric field is sharp as a 
function of frequency. This result is hardly surprising as 
there is no disorder in this model. However, it will serve 
as the basis for the perturbative treatment in the next 
subsection, and so is useful to analyze in some detail. 

To introduce uniform pinning in the model, one only 

needs to set hi = — J2r^o 2W ~ ^0 for all the sites i in 
Eq. |l^. The representation of the spin wave Hamiltonian 
in terms of phonon degrees of freedom, and its diagonal- 
ization, are formally identical to the steps used in the 
last subsection, provided one makes the replacement 

F{k)^ F{k)+vo = Fp{k). 

The resulting (pseudo)spin wave for this model is thus 

El ^F,(kY-A\G{kW. (24) 

In the limit fc ^ 0, EZ-^vq. Thus for collective modes 
in which all the electrons move together, only the center 
of mass degree of freedom is relevant: electron-electron 
interactions may be ignored, and one obtains the single 
electron excitation frequency |^l|. The long wavelength 
dispersion of Eq. |2^ may be obtained by evaluating F{k) 
and G(k) using an Ewald sum technique [Q. One finds 

F{k) « 2'Ke^llpok + 0{k^) 

G{k) w ne^llpo{k^ -f ikyf/k + O(k^), (25) 

(26) 

so that 

El^vo + 2Tre'^llpok + 0{k^). (27) 

The linear dispersion of the pinned collective mode spec- 
trum is purely a result of the long-range nature of the 
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Coulomb interaction. For a short-range interaction, the 
collective mode disperses instead as fc^ from vq. The lin- 
ear dispersion, and hence the long-range nature of the 
Coulomb interaction, turns out to play a crucial role in 
allowing the sharp response of the uniformly pinned sys- 
tem to survive the introduction of disorder. This will be 
discussed more carefully in the next subsection. 

Because one may exactly diagonalize the Hamiltonian 
in this model, it is convenient to compute the power ab- 
sorption using Fermi's Golden Rule. For an electric field 
of the form E{t) = Re Eo{x + iy)e*"*, one obtains for 
zero temperature 



Pico) 



Y.En\<n\ll 



(28) 



where |0 > is the groundstate of the system (no spin 
waves), and \n > represents the set of single pseudospin 
wave excitations (which in the present case may be la- 
beled by k rather than n) . For the uniformly pinned sys- 
tem, the matrix element entering Eq. ^ is nonvanishing 
only for = vq, so that one obtains a delta function 
response ai uj — vq, as expected. 

Eq. ^ is a good starting point for a perturbative treat- 
ment of disorder effects. In particular, the weight arising 
from the matrix element < 'T'|7j:^g|0 > must remain pro- 
portional to the system size in the thermodynamic limit 
for one particular mode n if the system is to retain the 
sharp response observed in Sec. || for an arbitrarily large 
system. In the next subsection, we show that, at least 
for weak disorder, this is indeed the case. 



we will see that Lc is extremely large compared to real 
sample dimensions. 

The quantity we will use for our perturbing parameter 
is the deviation of the pinning potential from perfect uni- 
formity, AUi, which is formally defined by the equation 
hi — — J2r^o 2iW ~ ^0 — AUi- Note from its definition 
that AUi = 0, a property we will use below. The 
power absorption (Eq. |2^) may be written via standard 
many-body manipulations |pT[| in terms of a Green's func- 
tion, in the form 

P{uj) = -e%2£;2^/m{^^d;i(*)dO(jrG;..(zj;c^-Hz<5)}. 



(29) 



The Green's function matrix entering above is given in 
imaginary time by 



G,A^j;T)=-<Trb\^\T)bf^\0)> 



(30) 



with ^.,v ~ 1,2, b^p — hi, and b^^ = 6-. The c-numbers 
dP^{i) entering Eq. 29 are those that diagonalize the 
pseudospin-wave Hamiltonian Eq. |2^ for the A: = 
mode in the uniformly pinned case (i.e., di{i) 



,k=0 



0, d2{i) = Ufc^o — Ij with w, v given by Eq. [2l|j. Us- 
ing the method described in Section |l| and Eq. ^0|, the 
equation of motion for the Green's function in imaginary 
time is found to be 



Gii{ij;iu!n) 
G2i(ij;«w„) 



C. Weak Disorder: Perturbative Treatment 



-1 
1 



Gii{ij;iu!n) 

G2liij;iU!n) 



In this subsection, we will formulate a perturbation 
theory for Eq. |2^, it terms of deviations of the pinning 
potential from uniformity. The point of this analysis is 
to understand how the sharpness of the response might 
survive the introduction of disorder. Towards this end, 
we will focus on the height of the 5-function response 
found in the last subsection, and develop an expression 
to the lowest non-trivial order in perturbation theory to 
see how much it is decreased by disorder. The result- 
ing expression, when disorder averaged, will turn out to 
have a formal, logarithmic divergence, which we interpret 
to mean that there is no true J-function response in the 
thermodynamic limit. However, the divergence is in fact 
cutoff by the system size, and we will see that even a 
conservative estimate of the integral for real system sizes 
indicates that it is in fact small, so that the perturbation 
theory is valid. The integral allows us to define a length 
Lc that is the characteristic size scale for electrons mov- 
ing together in phase in the lowest collective mode, and 



E 



where 



—Fp{Ri — Rk) —2G*{Ri — Rk) 
2G{Ri — Rk) Fp{Ri — Rk) 



^R.R' 



Gnikj;iujn) 
G2iikj;iuj„) 



(31) 



(1 



G{R) = 



4i?3 



(32) 



Note the the quantities in Eq. ^ are just the discrete 
Fourier transforms of the quantities Fp{k) and G{k) de- 
fined in the last subsection. The matrix elements G12 and 
G22 may be found by using Guiij; iw„) = Gli{ij; -zw„) 
and G22{ij]ii^n) = Gl^{ij; -iun)- 
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Eq. ^ may be solved in a manner closely analogous 
to that of Section ||: one needs to solve the eigenvalue 
equation 



E 



-AUAk-Fp{R,-Rk) -2G*{R,-Rk) 
2G{R, - Rk) AUAk + FpiR, - Rk) 



(33) 



and then the solution may be expressed in terms of the 
eigenvalues and eigenvectors. It is not difhcult to prove 
several properties of the solutions to Eq. p3 based on 



the symmetry of the matrix |39|. In particular, one may 
easily show that the eigenvalues come in pairs of equal 
magnitude and opposite sign, iwj. The eigenvector of 
the solution with negative ujj may be found from the 
solution with positive LOj with the transformation 



S^\k) = 4^''*{k), d)^'{k) = d\^>'*{k). 

Because the matrix being diagonalized is not Hcrmitian, 
these two eigenvectors are not in general orthogonal to 
one another. If one imposes the normalization condition 

J2ki\d2^ (^)\^ ~ \di\k)\'^) = 1 for positive eigenvalues 
{ujj > 0) and insists that the eigenvectors associated with 
each pair of eigenvalues of the same magnitude obey the 
above relation, the inverse of the eigenvector matrix may 



for u > 0. 

For a finite size system, Eq. ^ describes absorption by 
the system into a discrete set of states, with a weight that 
may be interpreted as the square overlap of the mode 
being excited with eigenvector of the fc = mode for 
the uniformly pinned case. In most situations, one ex- 
pects as the thermodynamic limit is approached that the 
weight associated with each mode vanishes with increas- 
ing system size, while the density of modes increases, to 
generate an absorption curve that is a continuous func- 
tion of frequency. What we have found in Sec. |l| seems 
to imply that the lowest mode in spectrum of the dis- 
ordered system retains a much larger overlap than all 
the other modes, even for large systems. To understand 
how this arises, we need to compute the overlap sum 
[j:^dU^)ei'\iy][j:k4ikrd^i\ky] in Eq. g for the 
lowest energy mode, and see how it scales with increasing 
system size. Although an exact calculation of this quan- 
tity is not possible for an arbitrary disorder strength and 
arbitrarily large system, we can at least estimate it for 
weak disorder to see if and when a finite overlap survives 
in the thermodynamic limit. 

Towards this end, we compute the eigenvectors of the 
matrix in Eq. ^ to second order in AUi . The method by 
which this is done is identical to standard non-degenerate 
perturbation theory in quantum mechanics [^0| . Denot- 

ing (4^")(i?i),4^")(i?2),...,4'^(i?i),4'\(^2),...) ^ VW, 
we may write the result of the calculation as 



be explicitly constructed. Writing this as e)i'{k) so that ■y(j) — \i_ jji '\\^/'^Y^~ 



E^,fceJf^(A:)d^/'^(fc) = <5,-,., wefind 



jU)* 



jU)* 



E 

k i^i) 



V, 



(fc) 



< k\AU\j > 



(0) 



OiAUf 



(35) 



for LOj > 0, and 



e(^")(fc)=4^>(fc), e^i\k) 



where V|,^'\ E^ are the eigenvectors and eigenvalues in 
the absence of the perturbation. The matrix element is 
given by 



for ujj < 0. With this explicit expression for the eigen- 
vector matrix inverse, it is possible to write the solution 
to Eq. 31 as 



< k\AU\j >= W, 



(fc)T 



^AU 
AU 



V, 



4'^Wei^'^(fc) 



for {^jLv) = (11) and (21). For = (22) and (12), 

the expression is the same as above, except one needs 
to take the complex conjugate of the numerator. Com- 
bining this with Eq. £91 and noting that d\{i) = so 
that E.fc Eu. dl{i)ef{iYd\^\k)*dl{k) is purely real, it 



with Wq'^''^ representing a column vector of the form 

{e^f{Ri),e^\R2),...,e^i\Ri),e^i\R2).-), and AU in 
the above matrix is an A'^ x iV diagonal matrix with en- 
tries AUi on the diagonal. The relevant quantity for our 
purpose is the "depletion " D{j), which one finds to be 



E 



< k\AU\j > 



E, 



(0) 



E 



(0) 



+ 0{AU^). 



(36) 



follows 



i i 

X [J2d°2{k)*d^f{k)*]6{Lu~LUj). 



(34) 



The zeroth order eigenvectors and eigenvalues are easily 
evaluated, as these were already essentially found in the 
last subsection. Thus the eigenvectors without disorder 
are conveniently labeled by a wavevector q and a sign ± 
denoting whether the positive or negative eigenvalue for 
a given wavevector is being referred to. Thus we use the 
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eigenvalues = ±E^ = ±[Fp{q)'^ - 4|G(g)|2]i/2^ and 
the corresponding eigenvectors have the form (i) 
-re" 



L-^^-^-t;,-, di'^ii) = -^e^«-«'u,-, with 



Fp{q) + Eg-- 



2E^ 



1/2 



mm 



/2E^ 



Fp{q) + Eg- 



1/2 



10,-/2 



(37) 



for E° > 0, and 




0,-/2 



(38) 



for Ej < 0. Since the factors appearing in the power 
absorption (Eq. are just the non-zero part of the 
eigenvectors above for g = 0, it is easy to verify that the 
weight associated with excitations into the fowest collec- 
tive mode in power absorption is just 1 — D{Q) (i..e, there 
will be a contribution to the sum over modes appearing 
in Eq. |3| proportional to 1 - -0(0).) Provided D{0) < 1, 
this will remain a (5-function contribution at zero tem- 
perature, even as absorption into the other modes may 
merge into a broad background. Thus, -D(O) measures the 
depletion of the sharp pinning mode response found for 
the uniformly pinned model. It should be noted that be- 
cause of our choice of system around which we are doing 
perturbation theory, there are no corrections to the pin- 
ning mode frequency at 0(AC/); the first non-vanishing 
correction is of 0(A[/^), so that the frequency of the pin- 
ning mode remains close to its uniformly pinned value. 

Using the forms for the uniformly pinned model in Eq. 
pq, one finds 



\u±\^\AU{q)f^ 



mm' 



(39) 



with AU{q) = AU.^e"'-^' . The second term in Eq. 

|39| is always finite and is small provided |A[/(g)/i?gp is 
small. The first may potentially diverge even for small 
|AL/((7)| because E^ — > uq as (f ^ 0, so that there is 
a vanishing energy denominator. However, for a given 
disorder realization, by our choice of the system around 
which we are performing perturbation theory, we have 
I AC/((j)| oc q' for small q. The energy denominator, using 
Eq. |2^, behaves as [vq — E^Y ~ (STre^/Qpo?)^, so that the 



integral remains finite, and can be small if AU{q) is small 
enough. 

This is the central result of this section, and several 
comments are in order. Firstly, the result of Eq. ^ 
is only finite in this analysis because of the long-range 
nature of the interaction. For short-range interactions, 
[vo — E^Y ^ and one ends up with a divergence 
no matter how small |A[/(g)| might be. Such a diver- 
gence indicates that one cannot stop at second order in 
the perturbation as we have done here, and that some 
self-consistent treatment is called for . Under these 
circumstances, one expects the resulting response to be 
broad, as is the case for most pinned CDW's As men- 
tioned in the Introduction, this is borne out by the QHA 
of Section ||, which shows that a broad response is indeed 
obtained if one uses a screened rather than long-range 
Coulomb potential (cf. Fig. 2). Secondly, in many calcu- 
lations where one averages over disorder configurations, 
the simplest choice of disorder models (white noise) intro- 
duces fiuctuations at all length scales, so that the disorder 
average \ AU{q)\'' = |A[/p is independent of wave vector. 
Thus, a disorder average of Eq. ^ would eliminate the 
zero in the numerator as (f— > 0, and introduce a logarith- 
mically divergent depletion in the thermodynamic limit, 
signaling a broadened rather than sharp response. 

In practice, however, we find that for physically rele- 
vant systems that the response remains sharp. The rea- 
son for this is that the divergence only arises for truly 
infinite systems. We can define a length scale Lc above 
which the depletion D{0) is of order 1, so that the S 
function response becomes significantly broadened. To 
do this, one must solve the equation 



1 

1 = — 



d'q 



PO A>27r/Le (27^)- 



(40) 



We can break up this integral into singular and non- 
singular parts as Lc ~^ oo, and so write Eq. ^ in the 
form 



1 = 



1 



2Tie 



dq \AU\^ 



27rpo 727r/Lc 9 



where we have used Eq. |2^. In the above equation, £ 
represents a length scale above which the collective mode 
dispersion for the uniformly pinned system is accurately 
represented by Eq. (i? « lOao would probably be 
sufficiently large), and ry(i?) represents the non-singular 
contribution to Eq. One may now solve for Lc, with 
the result 

f (27rpo)3(e;o)4 
Lc = ^exp| === [1 - ri{£)]} 



|AC/|2 



(41) 
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In Eq. |AJ7p is the disorder potential strength \ AU^\ 
written in units of c'^/kIo. For weak disorder, ri{£) <C 1. 
We wiU discuss the interpretation and consequences of 
Eq. ^ further in Section ^ below; for now, however, 
we point out that an estimate of Lc for physically rele- 
vant parameters shows that it is extremely large, much 
larger than the physical dimensions any real sample. This 
means that in practice, the depletion D{Q) will be small 
for weak disorder, since the sample size cuts off the di- 
vergence at length scales much smaller than the one at 
which broadening becomes significant. 

It is interesting to consider what would happen if one 
were to consider a model in which the displacement of 
electrons from their lattice sites was included as an effect 
of the disorder. The essential change is that the elec- 
tron centers Ri appearing in Eq. |3l| are no longer on 
lattice sites. For weak disorder, if one neglects lattice 
defects such as dislocations and disclinations, the effect 
of the lattice deformation may be described by writing 
Fp{Ri — Rj) — > Fp{Ri,Rj) — Fp{Ri — Rj) + SFp{Ri, Rj) 

&nd GiRi-Rj) G{R^,Rj) = G°{R^~Rj) + SG{Ri,Rj) 
in Eq. where Fp and are the couplings between 
electron lattice sites in some perfect reference lattice. 
One then may treat SFp and SG as perturbations in pre- 
cisely the same manner as we treated AUi above. The 
resulting depletion D{j) has precisely the same form as in 
Eq. with the only difference being that the matrix ap- 
pearing in the definition of {k\AU\j) now has off-diagonal 
matrix elements due to SFp and 6G. The resulting ex- 
pression for the depletion of the lowest collective mode 
D{0) has the same form as Eq. |9[ although the pre- 
cise form of |AC/((;)p will be more complicated than for 
the diagonal disorder model. The important point, how- 
ever, is that the energy denominators in the perturbation 
theory will be unaffected, so that one still expects only 
a weak logarithmic divergence if the disorder is not too 
strong, and a resulting sharp electromagnetic absorption 
for a finite size system. This demonstrates, albeit a pos- 
teriori, that our assumption of a diagonal disorder model 
does not alter the qualitative physics of the system, at 
least for weak disorder. 

As is clear from the above discussion, the precise re- 
sults for the length scale Lc depend upon the disorder 
model one uses for AU{q). In the next section, we intro- 
duce a simple model for this due to imperfections in the 
interface upon which the 2DEG resides. 

IV. INTERFACE PINNING MODEL 

In both the QHA and the perturbative analysis de- 
scribed above, one must make a specific choice for the 
(disorder) pinning potential. In the Section ^ below we 
present results for three different models. The first is 
a simple white-noise type model, in which some fixed 



fraction of the electrons are given an excitation energy 
vo + AUi = Up between the m = and m = 1 orbital 
states, and all other sites are unpinned. The simplicity of 
this model allows a clear comparison of results obtained 
from the QHA and the perturbative analysis. A second 
model we investigate is based on the idea that charged 
impurities are likely to be present in the spacer layers 
of heterojunction systems, and that some of these impu- 
rities, if close enough the the 2DEG, will substitute for 
electrons in the lattice [^. By choosing a small num- 
ber of lattice sites to have extremely large values of Up, 
we can model this type of disorder within the QHA. As 
will be seen, this leads to a band of collective excitations 
well above the phonon band, corresponding to localized 
excitations of the strongly pinned electrons. At low exci- 
tation frequencies, these electrons are essentially immo- 
bile, and so behave as if they were not degrees of freedom. 
This is by definition a strongly pinned system; yet we will 
see in Section ^ that the resulting pinning frequency is 
extremely small, and that an unacceptably large number 
of such charged impurities must be present in the system 
to account for the experimentally observed magnitude of 
the pinning frequency. 

We thus focus on a pinning mechanism which, to our 
knowledge, has not been previously discussed in the con- 
text of the magnetically induced WC: interface disor- 
der. The interfaces between GaAs and AlAs at which 
the 2DEG's reside are by design of very high quality in 
samples such as those of Refs. ||l^ and Nevertheless, 
they cannot be perfect, and it is generally accepted |2|] 
that such interfaces can only be defined to within a sin- 
gle lattice constant of the host semiconductors. A simple 
idealization of this is to model the interface as a series 
of pits and/or terraces, with the height of the interface 
fluctuating randomly up and down as illustrated in Fig. 
5. The typical scale of the pits and terraces formed by 
the imperfect interface are generally thought ||2^ to be 
of the order of several tens of Angstroms. For the hole 
samples of Refs. jl^] and the interface has an addi- 
tional corrugated structure with a size scale of 32A and 
a depth of IO.2A. 

There are several reasons for believing that this form 
of disorder may be important in the experiments of Refs. 
|]l5| and |l^. Firstly, the magnetic field at which one 
first sees a clear resonance at the lowest available tem- 
peratures is roughly 8T, corresponding to a magnetic 
length of ^0 ~ 90A, a very reasonable length scale below 
which one might expect the interface structure to trap 
electrons. If we model the pits in the interface as flat 
depressions of depth Az as in Fig. 5, one can estimate 
the trap potential for a large pit by assuming the electric 
field between the 2DEG and the charged remote donors 
is uniform |^]. In this case the the potential gained 
from placing an electron in a pit that is much larger than 
the magnetic length is AV — 2Trpoe^Az/K ~ A.35K in 
temperature units; k « 12 here is the dielectric constant 
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of the GaAs host for the electrons, and we have taken 
Az = lOA. It should be noted that for experimentally 
accessible magnetic fields, a typical pit size will gener- 
ally be much smaller than the magnetic length, so that 
the pinning energy of an electron trapped in a single pit 
will be of order ISVs^ /Iq, considerably smaller than the 
maximum possible value of A^. 

Because of the roughness of the interface, the WC in 
general will distort slightly so that some or all of the 
electrons may take advantage of the interface potential. 
Ultimately some fraction of the electrons will find equi- 
librium centers for their (groundstate m = 0) Gaussian 
orbitals that are particularly low in energy. Because the 
first excited {m — 1) state of each orbital is spatially far 
from the center of the m — Q state (in the sense that Iq 
is larger than the average pit size sq), the former state 
will not be correlated with the disorder potential, and 
thus is not on average lowered in energy as is the m = 
state. It should also be kept in mind that for a given 
groundstate configuration, the energy of an excitation of 
a single electron from the m — Q state to the m = 1 state, 
keeping all the other electrons fixed, will have its largest 
contribution from the electron-electron interaction rather 
than from the interface disorder. Provided the disorder 
is not too strong, the potential well in which an indi- 
vidual electron resides is thus to a first approximation 
circularly symmetric, so that the expansion in terms of 
angular momentum states used in this work is sensible. 

To fully specify the interface disorder model, we need 
to find the energy due to the disorder to excite an elec- 
tron out of its groundstate orbital (m = 0) into the lowest 
excited state (m ~ 1) for each site i, ei{m = 1) — ei(m = 
0) — vq + AUi. The distribution of AUi depends upon 
the assumptions one makes about the surface morphol- 
ogy of the interface, as well as whether other pinning 
sources besides interface roughness are present in the sys- 
tem. A fuller accounting of various possible models will 
be discussed elsewhere here we will focus on one 

reasonable possibility that reproduces both the magnetic 
field and density dependence of the resonance frequency 
observed in the experiments of Ref . . 

In this model we make the plausible assumption that 
the strongest pinning in the hole samples comes from de- 
fects and disorder in the interface corrugations of these 
systems. (It may be shown that perfectly regular corru- 
gations of this sort lead only to very weak pinning [^ .) 
As a simple model of this, we assume that the surface has 
pits of size scale sq = 30A and depth Az = lOA with a 
surface density n^. We will assume the average distance 
between pits 1 / is smaller than the nearest neighbor 
separation of the electrons a but larger than the magnetic 
length Iq. This means that each unit cell of the WC con- 
tains several pits in which an electron might be trapped, 
but when trapped an electron wavepacket covers one and 
only one pit. In what follows we will choose the density 
Ui to fit the experimental data at the lowest fields for 



which a resonance is observed; we then need to go back 
and check that the resulting pit density is consistent with 

To estimate the density of pinned electrons, we use the 
collective pinning approach that has been very success- 
ful in the context of vortices in superconductors and 
charge density waves [0,|l3|. The basic idea is to find 
a size scale Rc of correlated domains, whose positions 
may adjust more or less independently of one another to 
take advantage of the pinning potential, at the expense 
of lattice distortion energy. The average number of pins 
that may be found under electron wavepackets in a sin- 
gle domain if placed at a random location is given by 
Npin — PoR^T^l^ni] if we vary the position of the domain 
a distance of order the disorder correlation length (in 
this case, l/y^), one expects fluctuations in the num- 
ber of pinned electrons of order y/Npi^. Since the energy 
gained by an electron wavepacket when sitting over a 
pit is AVsq/Iq, in analogy with the Fukuyama-Lee-Rice 
model for charge density waves, the energy per electron 
gained from the disorder potential Up^^ is 



,FLR , 



AV 



IqRc Po 



1/2 



The introduction of such distortions inevitably leads to 
a lattice distortion of amplitude l/y/Ki over a distance 
Rc- Following Rcf. |l^, we assume the distortions are 
purely transverse, and so estimate the distortion energy 
per electron to be 



,FLR 



n^RlpQ 



where ^ is the shear modulus of the lattice. Minimizing 



the total energy 



FLR 



,FLR 



with respect to Rc, we 



arrive at the optimum domain size 



Rc 



2fllQ 



AVsUnnfpoy/^ 



(42) 



The resulting pinning energy per particle is found by sub- 
stituting this value of Rc into our expression for Up^^, 
leading to the expression 



,FLR 



To use this, we need an explicit form for the shear mod- 
ulus. This can be deduced from the transverse phonon 
dispersion relation of the WC in the absence of a mag- 
netic field, which has the form LOriq)^ — m*pa ' ^^i^^^ 
m* is the electron effective mass. The transverse mode 
eigenfrequency has been shown in the classical har- 
monic approximation to be given approximately by 

with Wp = A-Ke^ / \/?>m* no? . Using the above equation, 
one arrives at the estimate 
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/i w 0.3- 



(43) 



The resulting pinning energy per particle is finally given 

by 



,FLR , 



5.24- 



p272 



(44) 



In the perturbative regime, where the effective depth of 
the individual pinning centers, A1^Sq//q, is much smaller 
than the bandwidth of the phonon density of states, the 
pinning frequency is just given by the average binding 



u. 



FLR 



P 



with 



ut, given by Eq. Note within this approximation. 



energy per site, so that one finds ujpr 

FLR 

p 

i^pin oc Iq ^ oc B, which leads to a pinning frequency that 
increases with magnetic field, as seen in experiment. It 
is also interesting to note that ujpin oc oc which 
is also consistent with experimental results . We will 
obtain quantitative estimates for vq in the next section 
using the above expressions. 

When the pinning potential for an electron trapped 
by a pit approaches the bandwidth of the phonon den- 
sity of states, the perturbative estimate above becomes 
quantitatively inaccurate, and one needs to perform a 
QHA calculation to obtain reliable results for ujpin- To 
do this, we need the probability that a given electron 
will be trapped in a pit, which in this approach is given 

2 

by u^-'"^/ AV(^). We are then led to consider a distri- 
bution of pinning energies in which most electrons are 
unpinned, and a small fraction are pinned with an effec- 
tive potential vo + AUi — AVsq/Iq. The resulting model 
is then formally identical to the white noise model men- 
tioned above. 



V. RESULTS 

In this section, we describe in more detail the results of 
our study. As has been emphasized, the electromagnetic 
responses computed in the QHA in all the models we 
have studied are qualitatively the same: one finds a single 
sharp line dominating the absorption spectrum. This 
line is very robust in that its weight shows no discernible 
decrease with increasing system size, and for models in 
which the pinning potential of the individual electrons is 
small compared to the width of the phonon density of 
states (Fig. 4), the frequency of the resonance occurs at 
the average excitation energy per electron, as expected 
from the perturbative analysis of Section |l|. We begin 
by describing in more detail the results for disorder that 
is in this sense weak. 



A. Pinning by Weak Disorder 

A particularly useful model for comparison of the QHA 
and perturbative analysis is one in which a fixed fraction 
n-pin — Npin/N of randomly chosen electrons is pinned, 
with each such site assigned the same pinning potential 
vq + AUi = Up. According to the perturbative analysis, 



the pinning frequency Ljp 



lUp, regardless of the 



precise distribution of pinned sites. Figs. 1 and 2 show 
typical results for this model, with half the sites pinned, 
using Up = COIc^/kZo. Fig. 6 illustrates the values of 
^pin as computed in the QHA for 5 different disorder real- 
izations each at several different values of wq -I- AUi = Up, 
for a smaller density of pinned sites (10%), and AUi = 
for all other sites, along with the prediction of the pertur- 
bative analysis (dotted line). For small Up the agreement 
of the two approaches is quite good. Furthermore, there 
is almost no variation in Upin for a fixed value of Up when 
it is small as the precise realization of the pinned sites is 
changed. The interpretation of this observation is that, 
for small Up, LUpin is determined almost exclusively by 
the restoring force on the WC when only the center of 
mass coordinate is moved with the lattice itself undis- 
torted. Lattice distortions do introduce some deviation 
in the resonance frequency, generally pushing it below 
the value expected from the lowest order perturbation 
theory, in a way that does depend on the precise disor- 
der realization as may be seen in the figure. However, 
for any given disorder realization, a single mode always 
dominates the response; broadening is never observed for 
a fixed disorder configuration. 

We note that with higher order corrections to the per- 
turbation theory of Section III, more accurate predic- 



tions of 



can be made for a given disorder realiza- 



tion. For example, one of the configurations in Fig. 6 

= 2.32 X 10-"^. (The num- 



with Up — 0.003 has 



pin 



OJ. 



bers for both 

The lowest order perturbation theory predicts uj. 



pin 



and Up here are in units of e / kIq.) 



,(0) 



pin 



3 X 10 at this density of pinning sites, which 
is in error by approximately 35%. Using the methods 



of Section [H, the second order correction may be com- 



puted for this particular disorder realization, with the 
result 



(^pin ~ -I- LOp^^ = 2.19 X 10 which is within 
5.6% of the QHA result, a considerable improvement. 
(Note that the first order correction in AUi, "^pii, pre- 
cisely vanishes in our perturbative approach.) Another 
significant point, as stated above, is that the weight of 
this single mode shows remarkably little size dependence. 
For example, using the parameters relevant to Fig. 1, 
the power absorption per electron at the frequency of 
the sharp peak in the inset of Fig. 1 is proportional to 
Xxxii^pin)i^pin/N = 0.006209, where N = 1024 for this 
calculation. For precisely the same system parameters, 
but iV = 529, one finds Xxx{^pin)^pin/N — 0.006200, a 
slight decrease for the smaller system size. This decrease 



,(2) 
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is almost certainly related to the fact that the disorder re- 
aUzations in the two calculations are inevitably different, 
rather than to any systematic increase with increasing 
system size. Nevertheless, this result illustrates that one 
cannot discern a decreasing weight in the sharp response 
with increasing system size for the values of N one may 
handle in the QHA. 

The important question remains: to what extent, and 
under what conditions, does this response remain sharp 
for systems of experimentally relevant sizes? To address 
this question, we turn to the perturbative treatment of 
Section III. It was shown there that it is convenient to 



compute the power absorption by starting from a uni- 
formly pinned state - i.e., a pinning potential that is the 
same for all electrons in the system - and computing per- 
turbatively the corrections to this absorption spectrum 
due to the fact that the pinning potential is not truly 
uniform. For the uniformly pinned system, it is not at all 
surprising that energy absorption from a spatially homo- 
geneous (time-dependent) electric field is dominated by 
a single mode: the collective modes have a well defined 
wavevector fc, and only the A; = mode can couple to the 
electric field. In principle this k = mode is mixed in 
among all the modes when disorder is introduced. Thus, 
as discussed more carefully in Section III, we are led to 
ask whether there is a finite overlap between the fc = 
collective mode for the uniformly pinned case and the 
lowest frequency collective mode in a disordered system 
in the thermodynamic limit. 



In Section [II we showed that this overlap may be writ- 
D, with D given to second order in 



ten in the form 1 
perturbation theory by 



D 



1 

Po 



"2^ 



l4nAC/(g)| = 



K-l^|A[/(g1|- 



(45) 



with AC/((7) = -^Y.i^Uie'^-^'\ and MJi is the devia- 
tion of site i from the site-averaged value of the pinning 
potential, wq 
given in Eqs. |3^ and 

tion for the uniformly pinned system, Eq. |24|. For small 
g, as shown above (Eq. 



Explicit expressions for ut and are 
and is the dispersion rela- 



Vo + 2Tre^llpQq + 0{q^), 



so that D diverges in the thermodynamic limit, 
\AU{q)\^ 



if 



jAJjp > as g — > 0, as is typically the 
case for white noise potentials such as the one studied 
here. The meaning of this divergence is that essentially 
all the weight of the q — mode for the uniformly pinned 
system has been depleted by the disorder from the low- 
est energy collective mode, and is distributed among the 
other collective modes. However, in practice D is only 
divergent in the thermodynamic limit, since the integral 
in Eq. n3 has an infrared cutoff qmm — Stt/L, where L 



is the linear dimension of the system size. The depletion 
becomes significant (£> « 1) for system sizes L > Lc, 
which in Section III was found to be (Eq. |l|) 



£exp| — X= — I, 
^ \AU\^' 



(46) 



where |A[/p is the disorder averaged square potential 
evaluated in units of c^/kIq for g — *■ 0, and 2Tr/i is a 
wavevector below which the small q expansion for the 
uniformly pinned collective mode spectrum (Eq. [27[ ) be- 
comes reasonably accurate; presumably i ~ lOao, with 
the interelectron lattice spacing. 

For system sizes L <ti L^, the depletion remains small, 
so that a single sharp resonance should still be present. 
This means that one may interpret as a length scale 
for which the electrons in an infinite sample move to- 
gether coherently. To estimate the width of the result- 
ing resonance, we note that the fiuctuations in the pin- 
ning potential averaged over a length scale of Lc will 
^[^±^^] (AC72)^/^ so that the effective Q of the 



be 



resonance would be If Lc is very large, clearly 

this leads to an extremely narrow resonance. In fact, 
even a conservative estimate of Lc shows that it is larger 
than the physical dimensions of any real sample. Taking 
(AC72)^/^ « lOujpin, so that AC/ « lOt^p^K^o/e^ « 0.014, 
is the observed ||l5|,|l6| resonance frequency 



where 

of 1.25GHz, and we used a magnetic length ~ 81^ 
appropriate for a lOT magnetic field. (We expect that 
I AC/| is considerably smaller than this.) At this field, the 
filling factor for the experimental densities is v = 0.22, 
so that Lc ^ lO^'^flQ, far larger than any real sample. 
The conclusion then is that Lc^ L in Refs. [ p"5p^ , and 
the electrons at zero temperature react essentially as a 
single domain in response to the (spatially uniform) elec- 
tric field. Furthermore, the depletion D is in fact small 
in spite of the formal divergence when L — > oo at second 
order in perturbation theory; this accounts for the sharp- 
ness of the measured response at the lowest temperatures 

Finally, it must be emphasized that this effect arises 
purely due to the long-range nature of the Coulomb in- 
teraction. For short range potentials, — vq ^ q^ , and 
the divergence in Eq. ^ is much stronger. This obser- 
vation leads to another interpretation of the result: for 
short-range interactions, the phonon density of states for 
a uniformly pinned system at small frequencies is much 
larger than is the case for long-range interactions. (In 
the former case the phonon density of states jumps at 
the band edge, whereas in the latter it rises linearly from 
zero.) Thus, if one displaces the center of mass of the sys- 
tem (i.e., creates a g = excitation), for Coulomb inter- 
actions there are very few states into which the disorder 
can scatter this excitation. The calculations in this work 
demonstrate that the suppression of the phonon density 
of states at the band edge by the Coulomb interaction 
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is sufficiently strong to leave a finite oscillator strength 
of the pure center of mass (q=0) mode in the lowest fre- 
quency collective mode when disorder is included. 

The observation that Coulomb interactions are crucial 
to getting the sharp response is consistent with the re- 
sults of the QHA: Fig. 2 illustrates response for a sin- 
gle disorder realization when the Coulomb interaction 
is screened. For this calculation, we took an electron- 
electron interaction of the form 



v{q) = 2ne'^ / Ky/q'^ + ql, 

with qc = 2.QIq^, N — 225, and the parameters are oth- 
erwise the same as for Fig. 1. The broadening even for 
this relatively small system is already apparent. 



B. Pinning by Charged Impurities 

A commonly accepted model for understanding the 
magnitude of the depinning threshold observed in DC 
voltage-current measurements [p|, p^|j4^ ] is one in which 
charged impurities close to the 2DEG become incorpo- 
rated in the lattice as substitutions for electrons [ p6| . 
This is a paradigm for a strongly pinned system , where 
certain locations of the lattice are essentially "nailed 
down", and eliminated as degrees of freedom. In CDW 
systems, the pinning frequency for this type of disorder 
is estimated as the frequency of a phonon mode for the 
pure system, evaluated at wavelength q ~ 27r/d, where d 
is a typical separation between pinned sites. Using the 
QHA, we can quantitatively investigate this model to test 
whether it can account for the ~lGHz resonance as well. 

The method for emulating the charged impurity model 
within the QHA is to choose a small number of sites for 
which AJ7i is larger than the width of the phonon density 
of states. A typical collective mode density of states aris- 
ing from this type of disorder realization is illustrated in 
Fig. 7. Two peaks emerge, one at high frequency above 
the resonance frequency of the individual pinned sites, 
the other shifted slightly upward from the phonon den- 
sity of states for the pure system. It is easily checked that 
the motion of the electrons in the lowest energy modes 
leave the pinned electrons stationary, so that the magni- 
tude of the pinning potential on these sites is irrelevant 
for these modes; in practice these electrons are removed 
as degrees of freedom. 

Fig. 8 illustrates the frequency of the lowest collective 
mode as a function of pinned site density, with several 
different disorder realizations. As in all the results of 
this work, the electromagnetic absorption is dominated 
for each realization by only the lowest mode |47| . At low 
densities, the frequency is approximately linear in the 
pinning site density, a result quite reminiscent of what 
was found for the weak pinning model of Section VA 
above. The magnitude of the pinning frequencies found 



is surprisingly small; if one extrapolates the linear be- 
havior for small pinning densities, we find that ~ 13% of 
the sites would have to be pinned in order to achieve the 
1.25GHz resonance seen si v — 0.2. Such a high den- 
sity of charged impurity substitutions near the 2DEG is 
clearly inconsistent with the high mobilities these sam- 
ples exhibit in zero magnetic field. Furthermore, the pin- 
ning frequency is a monotonically decreasing function of 
field in this model, contrary to experimental observation. 
It thus seems quite clear that this model cannot be the 
primary pinning source relevant in electromagnetic ab- 
sorption. 

The behavior of the pinning mode at low impurity den- 
sities is surprising in its similarity to the results of the 
weak pinning model above, given that this is intrinsically 
a strong pinning mechanism. Indeed, the results of a 
strongly pinned CDW estimate grossly overestimate the 
pinning frequency at the densities we have studied. For 
example, if we assume a pinned site fraction = 0.05, 
then for filling factor v = 0.2 the strong pinning esti- 
mate yields an expected pinning frequency greater than 
7 X lQ~^e^ / kIq, while the frequency found in the QHA is 
roughly 5 x \Q~^e'^ / kIq. This great disparity between the 
strongly pinned CDW estimate and the result of a more 
realistic calculation we believe exemplifies the limits of 
using CDW theory to quantitatively analyze properties of 
the magnetically induced WC. Indeed, we believe the ten- 
dency for the system at low npi„ to behave so much like a 
weakly pinned WC may be understood if one thinks of a 
strongly pinned electron as a charged impurity to which 
a vacancy in the WC has become bound. One then may 
think of the strongly pinned WC as a weakly pinned WC 
with point defects, and it is likely that a perturbative 



analysis of this system as in Section [II would explain 



the sharpness of the response found from the QHA. If so, 
the strongly pinned WC for the purposes of ac response 
is equivalent to a weakly pinned, defective WC. 



C. Interface Pinning Model 

We now turn to the results of the interface pinning 
model defined in Section IV. As a first estimate, we use 



the perturbative result ujp 



u^^^, with w^^^ 



given by Eq. 44 . As discussed in Section fV , we set sq = 
30l, Ay = 4.3if « 90GHz, and take a = 480A, and 
K = 12, as is appropriate for the experiments of Rcf. [fsl. 
To estimate the density of pits n.i , we set vq 



- ,,FLR 



u 



p 



IGHz at S = 8T [[H, for which Iq « 90A, and use Eq. 
^ to solve for rii] the result is rii « 2.5 x lO^^cm"^. 
The resulting average distance between pits is l/y/rii w 
200A, which clearly falls into the range of validity of our 
model, Iq < l/y^ < a. We note with this estimate of 
Tii, the (Fukuyama-Lee-Rice) correlation length Rc with 
the use of Eqs. ^ and ^ is found to be Rc ~ 7.5a. 
Fig. 9 illustrates the result of the perturbative estimate 
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(solid line) as a function of magnetic field for this set of 
parameters. 

The actual potential for sites that are pinned turns out 
to be large enough that the perturbative approach over- 
estimates the value of uipm noticeably. We thus turn to 
the QHA to get a more accurate estimate of the pinning 
frequency. As discussed above, this requires us to find the 
what fraction of electrons is pinned; for the above param- 
eters, this works out to n„i„ = Npin/N = u'^^^/AV^ ~ 

10%. It is interesting to note that, since u^^^ oc l^"^ , the 
fraction of pinned electrons is independent of magnetic 
field. Each pinned electron has an excitation energy of 

^2 

Up — AV-j^ for these parameters. We can then proceed 
with the QHA precisely as in the calculations of Sec- 
tion VA above. Using the parameters adopted for the 



perturbative approach, we find aX B = 8T in the QHA 
^pin — 0.84GHz, slightly below the experimental value of 
IGHz. This indicates that we should raise our estimate of 
rii, since its value was chosen to match the experimental 
result for this particular magnetic field. Wc find that the 
pinning frequency uopin 

wlGHz in the QHA if the fraction 
of pinned sites Upin is raised to 11%; this can be achieved 



if we assume a pit density of = 3.0 x 10 



11 



The 



results of this calculation are illustrated for several val- 
ues of magnetic field in Fig. 9. As in Section VA, the 



perturbative result overestimates the pinning frequency 
by an increasingly large amount as the value of Up at the 
pinned sites increases; this leads to the sublinear growth 



of < 



'pin 



with magnetic field that is apparent in the QHA 



result. It is interesting to note that experimentally the 
variation of ojpin with B is indeed found to be sublinear 



VI. DISCUSSION 

The primary result of this work, that the zero temper- 
ature response of a pinned, magnetically induced two- 
dimensional WC to a spatially uniform, time-dependent 
electric field is sharp, to our knowledge is unanticipated 
in the literature. As we have seen, the key reason for 
this result is the long-range nature of the Coulomb po- 
tential. Early studies of a two-dimensional CDW system 
in a magnetic field did note that for weak disorder 
potentials a calculation of the (static) CDW domain size 
diverges when a l/r interaction is included. In Ref. [ p^ , 
it was argued that this divergence maybe removed be a 
careful treatment of the different energy scales for lon- 
gitudinal and transverse distortions of the crystal. This 
treatment computed a lineshape for electromagnetic ab- 
sorption by assuming the system may be thought of as 
independent Fukuyama-Lee-Rice (FLR) domains, with 
randomness in the pinning frequencies of the individual 
domains satisfying Awpm/wpm 1. By noting that there 
must be a length scale over which the variations in do- 



main pinning frequency effectively decouple the domains, 
an estimate for the expected lineshape was found. 

The model of Ref. ||l^ is in fact not very different than 
ours. One could think of the individual electrons in our 
calculations as FLR domains, and the model investigated 
would basically be the same. However, in that work it 
was assumed that the length scale over which Coulomb 
coupling may be ignored is the same as the FLR length, 
leading to a resonance with Q ^ 1 . This reasoning works 
well in the absence of a magnetic field, because each do- 
main has two possible polarizations for their motion in 
a collective mode. For the lowest frequency collective 
modes, the domains can execute a transverse motion that 
avoids long-range density fluctuations, which are very 
high in energy. Such long-range density fluctuations do 
appear when the motion of the domains is longitudinal, 
and so such modes contribute to the phonon density of 
states at high frequencies. For the transverse modes, the 
coupling of motion among the domains is weak, and the 
approximations of Ref. |l^ ] make sense. 

In the presence of a magnetic field, however, it is not 
possible to separate modes into transverse and longitu- 
dinal; these modes are inevitably mixed. The magnetic 
field causes the domains to move in a circular fashion, es- 
sentially circulating around their effective potential wells. 
However, if different domains circulate at different fre- 
quencies, there will necessarily be long-range density fluc- 
tuations, so that such modes will be high in frequency. 
Low-frequency modes can be achieved if the coupling be- 
tween domains is explicitly included, so that correlations 
in the motions of different domains may be introduced. 
This means in a magnetic field the system will have a 
dynamical correlation length that is different than the 
static (FLR) correlation length, and in this work we have 
seen that this dynamical correlation length is extremely 
large. We emphasize that this behavior is very different 
than what occurs in most CDW systems, where the static 
and dynamic correlation lengths are basically the same 
||l3|,^. The new physics arises because of the unique 
combination of two-dimensions, the long-range Coulomb 
interaction, and a magnetic field, and leads to the sharp 
resonance found in this work. 

Experimentally, it is somewhat surprising that only the 
most recent measurements have uncovered the sharpness 
of this resonance, whereas there have been a number of 
earlier measurements of rf, surface acoustic wave, and 
microwave response in this system |p|- jl0| , p^ which found 
broad resonances. There may be several reasons for this. 
Firstly, in order to couple to the 2DEG, the experimen- 
tal methods probe the system at a finite k, not with a 
purely spatially uniform electric field. From a practical 
viewpoint this is necessary, as the oscillator strength for 
transitions at small k is very small, making detection of 
the absorbed energy difficult. However, a repetition of 
the analysis of Section [II shows that the absorption is 
sharp only for fc = 0: the second-order correction in per- 
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turbation theory diverges more strongly at finite k than 
for A: = 0, indicating that the response broadens as k 
increases. This observation is supported by the resuhs 
of Ref. for which several absorption peaks are ob- 
served, presumably representing harmonics of the funda- 
mental probing wavelength. The widths of these peaks 
are found to increase with increasing frequency. Thus, it 
may be that some of the methods used in previous work 
coupled to or mixed in large enough values of k to wash 
out the resonance. 

A second important aspect of the result is that it ap- 
plies only at zero temperature; clearly one expects ther- 
mal broadening of the resonance, which we have not ad- 
dressed in this work. This means that one can only 
expect such sharp resonances at the lowest tempera- 
tures. One study of the temperature dependence of the 
resonance indicates that it is only apparent below 
lOOmK, and that it sharpens very rapidly in the range 
50mK 30mK. It seems quite reasonable that thermal 
broadening is responsible for the absence of this sharp 
resonance in previous experiments, especially since the 
temperature below which the sharp resonance settles in 
may be sample-specific. 

A number of open questions remain unresolved by this 
work. Firstly, because this is a zero-temperature study, 
we have not been able to understand the detailed absorp- 
tion lineshape of Ref. jljjl^. Beyond thermal effects, the 
large scale coherence of the electron motion at zero tem- 
perature indicates that specific dissipation mechanisms ~ 
in particular, edge states - of the WC may prove 

important in this context. Beyond the lineshape, some 
aspects of the magnetic field dependence of the resonance 
remain unexplained: at the highest magnetic fields, the 
resonance shows little or no field-dependence, while the 
Q of the resonance continues to increase with field. This 
work gives a natural explanation for how the frequency 
may be an increasing function of field in an interface pin- 
ning model; however, it is not clear how one could obtain 
a field- independent resonance [?] . Investigations of these 
issues are currently being pursued. 

VII. SUMMARY 

In this work, we studied the response of a two- 
dimensional Wigner crystal in a strong magnetic field 
to a spatially uniform, time-dependent electric field, at 
zero temperature, pinned by a disorder potential. A new 
approach to computing the response functions of a local- 
ized electron system in the lowest Landau level was in- 
troduced, the quantum harmonic approximation. It was 
found that the response is sharp; i.e., there is a resonance 
that is not disorder broadened. For weak disorder, this 
effect was shown within perturbation theory to survive 
for macroscopically large samples, because of the emer- 
gence of an extremely large length scale Lc that repre- 



sents the distance over which electrons oscillate together 
in the lowest excited state of the system. The fact that 
Lc is so very large was shown to result primarily from the 
long-range nature of the Coulomb interaction. A model 
of interface pinning was shown to reproduce both the 
magnitude and some aspects of the field dependence of 
the resonance as observed in experiment. 
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APPENDIX A: COMPUTATION OF 
INTERACTION MATRIX ELEMENTS 

The interaction matrix elements t™jm2m3m4 '^^^ be 
greatly simplified by the use of the strict translational 
periodicity assumed in this work. In general, they may 
be written in the form 

X</'mi(^- Ri)<j)m2{r- Ri)<j)*n3{f~ Rj)(j),ni{r ~ R-j), 

where (j)mir) is the mth angular momentum state cen- 
tered around the origin. The quantity v{q) is the Fourier 
transform of the electron-electron interaction, which for 
most of this work will take the form 27re^/g. (To de- 
scribe a screened Coulomb potential, one may take v{q) = 
27Te'^ / y/q^ + q^. This was the form used in the calcula- 
tions leading to Fig. 2.) Because we are imposing pe- 
riodic boundary conditions, the interaction needs to be 
replaced with one that is periodic in a superlattice of 
unit cells, each with N electrons, whose positions inside 
each cell are identical. This is accomplished by making 
the replacements q ^ G, J d^q — > the above 
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expression, where {G} are the reciprocal lattice vectors 
of the superlattice, and Vc is the supercell area. Making 
this replacement, and switching over to bra-ket notation, 
we have 



^ mim2'm,3m4 



(Al) 

Although for the unscreened Coulomb interaction v{G) 
diverges for G = 0, we may formally set it to zero if 
one includes the effects of a uniform neutralizing back- 
ground. The matrix elements < mi\e~'^'^'^\m2 > may be 
computed with some work analytically. The result is 



example, our calculations with N = 1024 electrons in- 
clude approximately 60,000 different values of G in the 
sum. Because we include so many of these, the QHA can 
accurately reflect what happens at length scales shorter 
than a magnetic length, allowing one to treat disorder 
potentials that vary on a short length scale. This is the 
advantage the QHA has over the method of Ref. [ p4[ , 
which is typically limited to several hundred values of G. 
In the text, we took advantage of the sum rule 



"^mim200 — / . "^miT 



^mim200 ^ Smim2 (A-3) 



which we now demonstrate. We begin with simple obser- 
vation that 



< mi|e ^^"^\m2 >= 



rmi-m2 



7712! 



7771 ! 



(G/o)2 



(mi-?n2)/2 



where the {g} are the reciprocal lattice vectors of the 
electron lattice (i.e., not the superlattice). Then 



for 7771 > '772, where 6g is the angle between the the 
vector G and the x axis, and L'ml~"^^ is an associ- 
ated Laguerre polynomial. The expression for 7772 > 
mi may be obtained using < 77ii|e~*'^ ''|7772 >=< 
7772|e-*(-'5)-f|j^^ >* 

suit in the form 



It is convenient to write this re- 



< 777i|e 



7712 >= e 



i(T?i2 — rni)(eg + .j) 



2> p 



.(G), 



From Eq. A2 it is clear that < 0|e''^ ''|0 > depends 



only on the magnitude of g and not its orientation, so 
that the orientation angle enters the sum above through 

Since each reciprocal 



< 7771 e" 



-ig-r 



7^2 >oc e 



i(m2—mi)0g 



lattice vector has five others of the same magnitude ori- 
ented at angles that are integral multiples of 7r/3 away 
from 9g, if follows that the phase factor will cause the 
above sum to vanish unless mi ~ 7772. Thus 



with 



.(G) = 



7772! 

7771 ! 



(Glof 



2 -1 (mi-m2)/2 



^ ] ^mim200 ^ <5mim2- 



(A4) 



Similarly, for [/, 



mim2 00 



we have 



-Gll/2 



for 777i > 7772. For 777i < 7)72, the cxprcssiou for i^mj^mj (G) 
has the same form as above, only with the indices 7711, 7772 
interchanged. 

The interaction matrix element can now be written in 
the form 



^ 7nim2m^rn4 



U: 



mim200 



7;(G) < 7ni|e-'<= nr772 >< 0|e''- ''|0 > 



iG-ri 



G 



As above, the quantity entering the sum involves the ori- 
entation of G only through e«(™2 -"11)6(3^ gQ ^}jjg g^j^ must 



also vanish unless ?77i = 7712. Together with Eq. A4, this 



= - Vt;(G)e-^^-(«'-«^)^^„,„,(G)i^™3„,(G) 



proves the sum rule Eq. A3 



^ ^_2^ym4-m3^l"il-"i2| + |m3-™4|gi["i2-mi+m4-m3]8(5 

which is particularly convenient for calculations. For a 
given value of Ri — Rj one in practice can include a very 
large number of reciprocal lattice vectors in the sum; for 
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FIG. 1. Center of mass response function Xcm for ^ system oi N = 1024 electrons a,t f = 0.2. Half the sites are chosen at 
random to have a pinning potential AUi = O.Ole^ /kIq (see text); the other half are unpinned. Two states per site were retained 
in this calculation (see text). Inset: Relative power absorption for the same system. 
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FIG. 2. Center of mass response function Xcm for a system of size N = 225 electrons at z/ = 0.2 with a screened Coulomb 

potential of the form v{q) = l-Ke^ jnic^ + c^)^^"^ , with qc — 2.0/io. Half the sites are chosen at random to have a pinning 
potential tJJi — O.Ole^/K^o (sec text); the other half are unpinned. Two states per site were retained in this calculation. The 
result indicates that screening broadens the response, so that the sharpness of the response illustrated in Fig. 1 is related to 
the long-range nature of the Coulomb interaction. 
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FIG. 3. Collective mode density of states for an unpinned WC as calculated using the quantum harmonic approximation. 

N = 225, u = 0.2, and 5 states per site are retained in the calculation. The result illustrates that including higher order angular 
momentum states introduces collective modes at high energy, which may be interpreted as an analog of Wannier excitons (see 
text). Such collective modes have little effect on the low energy dynamics of the system, so that one may drop the high angular 
momentum states without introducing serious errors in the pinning response of the system. 
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FIG. 4. Collective mode density of states for an unpinned WC as calculated using the quantum harmonic approximation 
(QHA). A'^ — 1024, — 0.2, and 2 states per site are retained in the calculation. This has a form consistent with the phonon 
density of states for a classical, unpinned WC in a strong magnetic field, and demonstrates that the QHA produces sensible 
results for the test case of a WC in the absence of disorder. 
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FIG. 5. Schematic representation of the interface disorder model. The interface has pits and terraces so that the setback 
from the ionized donor plane varies by a distance of order Az. Electrons localized by the magnetic field into wavepackets in 
the X ~y plane of size scale ~ Iq may lower their energy by moving their centers close to pit centers, since in the area of the pit 
the electron will be closer to the donor layer. For larger pits or decreasing magnetic length, more of the electron wavepacket 
overlaps with the pit region, decreasing the energy of the wavepacket. 






+ 



+ + + + 



+ + + + + 



27 



FIG. 6. Pinning frequency LOpin for a system of N = 529 electrons, at filling fraction u = 0.2, as a function of the pit 

potential. 10% of the electrons arc pinned at random sites with energy AU = Up; the remaining electrons are unpinned. Stars 
are the results of the QHA for five different disorder realizations at each value of Up, with two states per site retained in each 
calculation. Dotted line is the expected pinning frequency as calculated in lowest order perturbation theory. For small values 
of Up the approaches agree quite well; for larger values the perturbation theory sets an upper bound on the pinning frequency. 
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FIG. 7. Density of states for a strongly pinned WC, computed using the QHA, with N = 529, v = 0.2, and two states 
retained per site. 10% of sites have a pinning potential vq + At/^ ~ Up = O-OAe^ / nla, a relatively large pinning potential; the 
rest are unpinned. A high energy set of collective modes may clearly be seen separated above the main peak. These modes 
arise due to localized collective modes of the strongly pinned electrons. For the modes in the lower, main peak, the pinned 
electrons are essentially stationary. 
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FIG. 8. Pinning frequency Wpi„ for a system oi N = 529 electrons with strong pinning centers, uipin was computed 

using the QHA with two states per site at fiUing fraction v = 0.2 as a function of Upin, the fraction of pinned sites, with 
vo + AUi = Up — OMe^ /kIo the potential of the pinned sites. (1 — Upin of the sites are unpinned.) Stars are the results for five 
different disorder realizations at each value of npi„. The computed frequencies are significantly lower than what is expected 
from an estimate based on strong pinning of an elastic medium. 
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FIG. 9. Pinning frequency ujpin in the interface pinning model as a function of magnetic field, for AV = 90GHz corresponding 
to a pit depth of Az — 10 A. The average pit size so = 30A, and the electron density in this figure is fixed at po = 5 x lO^^cm"^. 
Solid line illustrates the perturbative estimate of cupin for a pit density of Ui = 2.5 x 10^^ cm^^. The magnitude of the pinning 
frequency may be seen to increase linearly with magnetic field, due to the decreasing magnetic length of the Gaussian orbitals 
in which the electrons reside in the groundstate. Stars represent the results of a quantum harmonic approximation (QHA) 
calculation with the same parameters except for a slightly increased pit density, m = 3.0 x 10^^ cm~^. Results for five different 
disorder realizations are shown for B = 7, 9, 11, 13 and IST. 
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